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SI.  Intro  met  ion. 

In  nan.y  important  applications  of  mathematical  physics,  the  final 
numerical  answer  depends  upon  the  evaluation  of  a  Laplace  transform 

-00 

(1)  g(y)  ■  /  e_Xyf(x)dx,  Re(y)  >  0. 

Uq 


Integrating  by  parts  n  times,  the  result  is 


U) 


g(y) 


ILkl  .  il£l . 

y  2 

*  V 


♦ 


“  /’a0«-x=rf(n)(x)dx. 

y  Jo 


If  f(0)  /  0,  we  nay  use  the  first  n  terms  as  an  approximation  to  g(y). 
The  error  term  is  trie  integral  on  the  right-han i  'side  of  (,.)  which  is 

bounded  by 


(3)  y“(  •*!)  |f(n)(x)| 

CKxvcd 

for  y  real  md  by 

U)  |y|"n  Max  |f^'(x)l  /  Re(y)  , 

C<X<0D 

in  case  y  is  complex.  In  nar^y  cases  of  importance,  the  series 

s(y)  -  f  f(k)(0)/yk41 

k-0 


(5) 
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is  divergent  for  all  y.  Consequently,  in  order  to  make  maximum  use  of 
(2),  n  must  be  chosen  as  a  function  of  y  to  minimize  the  expressions 
in  (3)  or  (4).  The  fact  that  S(y)  is  a  divergent  asymptotic  series  as 
j —> oo  imposes  a  lower  bound  on  the  magnitude  of  the  error  term. 

The  question  arises  then  as  to  whetner  it  is  possible  to  obtain 
a  better  approximation  to  g(y)  than  that  given  by  (2)  using  some 
rational  function  of  l/y ,  with  coefficients  determined  by  the  f'  ‘'(0). 
One  way  of  obtaining  those  rational  approximations  is  that  of  expanding 
S(y)  into  a  formal  continued  fraction  and  using  the  n-t.h  convergents. 

The  purpose  of  the  present  note  is  to  indicate  another  method 
which  under  all  circumstances  yields  m  error  tern  leas  than  or  equal 
to  that  given  by  (3)  or  (4),  and  in  general  may  be  expected  to  improve 
it  considerably.  After  indicating  the  general  method,  we  apply  it  to 
the  two  most  commonly  encountered  non-elenentary  transcendents 


The  following  result  is  obtained: 

Theorem:  If  y  >  0,  there  exist  for  each  n  >  1,  polynomials  in  p,  Qn(y), 
R(y)  with  the  properties  that 


/0°£*y 

o'o 


lx 


R"(/)  nl  nl 

0^7  S  2*1.1  yn.l  Q  (r)  221.1  yn*l  * 


(7) 
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where 


(8)  (a)  Q^y)  =»  1  ♦  ^  ^y"1  , 


,n-k*l  n-k*l 


(b)  a.  -  nl  — — - j  -  . 

K  [(n  ♦  1  -  k)  l]*  /=1  \ 


“A^JL  If  •)  \ 

TT  ■ 


(<•)  Rn(y)  -  Pn(y)  ♦  Z  v'lpn-i(j,)  -  Vi*""'1  • 

i  *  1 


k-1 

T- 

ife 


(<*)  rk(y)  *  (-1)*  y”^*1^.  k  *  l,  2, 


If  y  is  complex.  the  ri^ht-sin e  of  (7)  is  replace^  by 


(9) 


nl 


^  !^)n  l^(y)|He(y) 


A  similar  result  holds  for  any  integral  of  the  form 


(10) 


e-Xy 

(1  ♦  x)a 


dx  • 


Applied  to  the  integral  in  (7),  I^(y),  the  method  discussed  in  (2)  et  seq. 

/  n*l 

yielss  the  sound  nl  /  y-  .  The  best  possible  bound  is  obtained  by  taking 
n  *  [yj»  yielding  an  error  tem  approximately  equal  to  e‘y//2*y  .  The  best 
possible  bound  obtained  from  the  above  theorem  is  obtained  by  taking 
n  ■  2  y  and  is  approximately  e”  ^  \A-ny  . 


Although  these  functions  in  the  cases  a  ■  l/2  and  1  have  been 
extensively  tabulated  for  real  y,  there  may  3till  be  some  practical 
application  of  the  inequalities  we  obtain,  since  in  several  important 
applications  theso  functions  occur  with  complex  argument.  We  discuss 
below  a  device  which  may  be  employed  if  Re(y)  is  small. 


§2.  The  method . 

Let  us  consider  first  vhe  simplest  case  where  f(0)  /  0,  y  is 

real,  and  we  will  be  satisfied  with  an  approximation  with  an  error  of 
2 

order  1/  y  as  y— >oo.  Integrating  by  parts  once,  we  obtain 
(1)  g(y)  -  ^21  ♦  ^  I  e-xy  IX(x)dx 

7  yd  0 

■  *  y  f0  ®"Xy[fl(x)  *  af(x)  * 

-f  P00 e"*7  f(x)dx  -  . 

7  dO  y 


From  tnis  follows 


U) 


whence 


(3) 


<  Max  |f^(x)  ♦  af(x)  ♦  b|  /y2  . 

(XxvOO 


The  parameters  a  and  b  are  now  to  be  chosen  so  as  to  minimize  the 
right-hand  side.  In  general,  this  is  a  difficult  problem,  which  we  do  not 
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propoae  to  discuss  here.  Nevertheless,  in  many  cases,  a  few  trials  will 
yield  values  of  a  and  b  which  appreciably  reduce  the  error  term. 

Let  us  now  proceed  to  the  general  case.  After  n-integration."  by 
parts,  we  have 


U)  g(y)  -  ?_(/)  ♦  ~  P  e"*7  f^(x)dx, 

n  yn  ^0 

where  we  have  set 


(5) 


(6) 


n-1 


Pn(y)  -  t  f(k)(0)y-(k,1)  , 

*•  I  .  .  A 


n  »  1,  2,  •• *. 


k-0 


We  now  write  g(y)  as  follows: 


g(y)  •  ?n(y)  ♦  J  e"*y^(n^(x)4  a^^^x)  ♦  •••♦  anf(x) 


♦  a 


n*l 


Utilizing  (u)  to  eliminate 
to  the  inequality 


(7) 


g(y)  - 


R„(y) 

0^77 


00  /  •  \ 

the  integrals  /  e""*^  f'n”  '(x)dx,  this  leads 

,J0 

Max  |f^  (x)  ♦  a^n"^(x) ♦  • • •  ♦  anf(x)  ♦  a  jJ 
0<x<oo 


< 


where  we  have  3et 


Qn(y)  ■  1  ♦  I  a  y”1 
i-1  1 

Rn(y)  ■  ?n(y)  *  t  ^  y'1  ?n_i(y)  -  Vi  *' 


(8) 
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The  error  term  depends  upon  the  polynomial  Q^(y)  and  the  functional 

(9)  cn(f)  -  Min  Max  |f(n)(x)  ♦  a,  f(n*1J(x)  ♦  •  •  •  ♦  af(x)  ♦  an  .  I 

OCX<00 

We  may  conceive  of  cn(f)  as  a  measure  of  the  deviation  of  f(x)  from  a 
function  of  the  form 

N  V 

(10)  h(x)  -  2!  Pk(x)«  k  » 

k-1  K 

where  each  p^(x)  1*  a  polynomial  in  x. 

In  the  succeeding  section  we  determine  cn(f)  for  the  case  where 
f(x)  -  l/l*x  . 


§  3 •  Application  to  the  Exponential  Integral. 

Referring  to  (6)  of  §1,  we  see  that  it  is  sufficient  to  treat  the 
Laplace  transform 


(1) 


Consequently,  we  wish  to  determine 


(2) 


Min  Max 
0<x<oo 


n); _  *],(**-  *)» 

(1  *  x)n 


( 1  ♦  x) 


n^l 


♦  •••  ♦ 


(-l)n  i 
1  ♦  x 


(-l)n 


n*l 


1  1  ♦  x  * 

Upon  setting  ^ ^ —  »  this  reduces  to 


(3) 


Min  Max 
-1<x<1 


nl  (l*x)n41  (n- 1)1  a1(l4  x)n 


-,n*l 


♦  •  •  •  ♦ 


(-l)n  a 


n*l 
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3ince  the  degree  of  the  polynomial  ia  (n*  1)  and  the  coefficient  of  the 

1 

highest  degree  tom  ia  nl  /  2  ,  it  follows  that  the  minimal  polynomial 

ia  the  (n+l)st  Chebychev  polynomial 


U) 


nl  *„.!<«> 
„n*l 


nl  coa(n  ♦  1)  arc  co3  x 


,2n*l 


.  ,,  vn+1  (n-  1)  1  a,  (1  ♦  x) 
nl  (1*  x)  '  1'  '  ,  ,xn 

*■ - i -  - - -  •••  ♦  (-1) 

„n*l  x  ' 


n+1 


and  that 


(5)  c  -  nl  /  2 

'  n  .1  ♦  x  /  ' 


2n*l 


Referring  to  (7)  of  the  previous  section,  we  see  that  before  we  can 
apply  the  inequality  with  confidence,  we  must  know  something  about  the 
location  of  the  roots  of  Qn(y)  ■  0,  Using  (4),  we  shall  determine  d,he  a^ 
and  show  that  they  are  all  positive,  whence  Q^Cy)  >  1  for  y  >  0. 

For  n  •  1,  2,  3,  the  roots  of  ^n(y)  *  0  are  all  negative,  and  we 
hazard  a  gues3  that  this  is  true  in  general.  This  result,  if  valid,  would 
be  of  importance  in  applying  the  inequality  to  the  case  where  y  is  complex. 

Since 


(6) 


W-1  *  (1*x»  ■ 


n+1 

I 

k-0 


‘k! 


(-1)  (1  ♦  x) 


kl 


we  obtain  from  the  identity  of  (4), 

(7) 


/  , vn-k*l  ,  , 

Vk*l  kI 


,n-k*l 


nl  *n+l 
n*l 


(k) 

T  (-1) 


kl 
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To  determine  T  ^  ( —1 ) ,  we  use  the  fact  that  Tn+^(x)  satisfies  the 
differential  equation 

(8)  (l-x2)T^'4l(x)  -  xT^x(x)  ♦  (n  ♦  l)2  Tn<>1(x)  -  0. 

Hence  T^|(x)  satisfies  the  equation 

(9)  (l-x2)T^*2)(x)  -  (2k*  l)xT^1J(x)  ♦  ((n*l)2  -  k2)r^j(x)  -  0. 


Thus 


(10) 


T&11  (-D 


(n*  l)2  -  k2 


2k  ♦  1 


<-»  - 


whence  finally 


(ID 


Ti*i  (-1) 


(-D 


n*k*l  -rr  ( (n  ♦  l)2  -  k2 


TT 

k»l 


2k  ♦  1 


k  -  1,  2, 


(-Dn+1  , 


k  -  0. 


Substituting  in  (7),  we  obtain 

,n*l-k  n*l-k 


(12) 


*k 


JSLL 


[(*»  ♦  1  -  k)i]‘ 


lf’k  f  (n  <1)2  •  /  \ 

2 1  l,  "  ♦ 1  >  ’ 


k  4  n  *  1  , 


nt  , 


k  -  n  *  1 . 


This  completes  the  proof  of  the  theorem  stated  above. 

If  Re(y)  is  small,  the  tern  1/Re(y)  may  greatly  increase  the  error 
term,  particularly  for  small  n.  'We  nay  overcome  tills  to  some  extent  by 
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writing 


,  00  /N  ,  ,00  -xy  [,  ,,n  ,  a,  (n-l)l  , 

(13)  «  '(x)dx  - --  /  - - -z  A - :r?*  •••♦(-!)  » 


y  CO  (l*x)  (l*x)  (l*x) 


minus  the  appropriate  terms. 


From  this  we  obtain  as  a  bound 


Min  Max  — ^y  ♦  •••  +  (-l)n-^  a  ,  /iy|n  -  -y  — r — —  , 

at  0oc<oo  (Lx)"*1  n~l  22n*1|y|n 


which,  if  Re(y)  is  snail,  will  be  less  than  nl/22n*hyln  Re(y), 


Application  to  the  Error  Function. 

Referring  to  (6)  of  6,1,  it  is  sufficient  to  consider 


s(y) 


r00  o'** 

- 172 

to  (l*x)1/id 


In  this  :ase,  we  have  the  problem  of  determining 


Min  Max  lf^(x)  ♦  a^f^n”^(x)  ♦  •••  ♦  a^f(x)  *  j  » 


0<x<ao 


where  f(x) 


Wm 

(l*x)  Taking  an>1  *  0,  we  obtain  as  a  bound 

-1/2  (^)  1^2^) 

VI  WW  U-w  <  1  I  1 _ _ . . 


Min  Max  (1  ♦  x) 

ai  O^X^CD 


(1  -  X)n 


♦  •  •  •  ♦  a 


<  Min  Mac  •••  t*»-l i 


a,  Oocwoo  2n  (1  ♦  x)n 


:  -175 

- 
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Similarly, 

“®  can  obt  dn  bounds  for  the  general  Integral 

(4) 

1  (y;  -  /  - -  dx  , 

L' 0  (lax)* 

or  integral 

s  of  the  fora 

» 

(5) 

^  n 

°o  -xy-x 

Ja(y)  '  /  ,  a  &  • 

■  0  (1-  X)a 

